% 2-D 4 group KOEBERG test problem, 
% ref keff =  1.007962, i get 1.007997 w/ extended 

it = 1;

dat = [ ...
    2.491869  0.067927  0.008228  0.745248  0.000000  0.000000  0.000000  0.000000
    1.045224  0.066508  0.000536  0.254328  0.063789  0.000000  0.000000  0.000000
    0.677407  0.070757  0.007058  0.000424  0.000486  0.064381  0.000000  0.001245
    0.375191  0.069235  0.083930  0.000000  0.000000  0.000003  0.050849  0.000000
    2.492653  0.067275  0.008295  0.745248  0.000000  0.000000  0.000000  0.000000
    1.049844  0.065296  0.000713  0.254328  0.063112  0.000000  0.000000  0.000000
    0.676610  0.070432  0.009230  0.000424  0.000478  0.063078  0.000000  0.001543
    0.379481  0.086595  0.108244  0.000000  0.000000  0.000003  0.048420  0.000000
    2.491978  0.066922  0.008285  0.745248  0.000000  0.000000  0.000000  0.000000
    1.051910  0.064628  0.000713  0.254328  0.062765  0.000000  0.000000  0.000000
    0.677084  0.069952  0.009214  0.000424  0.000473  0.062404  0.000000  0.001598
    0.381453  0.089675  0.108087  0.000000  0.000000  0.000003  0.047549  0.000000
    2.492535  0.067963  0.008459  0.745248  0.000000  0.000000  0.000000  0.000000
    1.045298  0.066632  0.000923  0.254328  0.063737  0.000000  0.000000  0.000000
    0.674684  0.072139  0.011714  0.000424  0.000486  0.064330  0.000000  0.001630
    0.374240  0.092630  0.133600  0.000000  0.000000  0.000003  0.049518  0.000000
    2.492329  0.066940  0.008409  0.745248  0.000000  0.000000  0.000000  0.000000
    1.051953  0.064694  0.000921  0.254328  0.062737  0.000000  0.000000  0.000000
    0.675683  0.070681  0.011675  0.000424  0.000473  0.062376  0.000000  0.001797
    0.380606  0.102043  0.134282  0.000000  0.000000  0.000003  0.046859  0.000000
    2.491521  0.066584  0.008400  0.745248  0.000000  0.000000  0.000000  0.000000
    1.054029  0.064020  0.000921  0.254328  0.062386  0.000000  0.000000  0.000000
    0.676197  0.070201  0.011651  0.000424  0.000468  0.061696  0.000000  0.001852
    0.382438  0.105135  0.133974  0.000000  0.000000  0.000003  0.046005  0.000000
    2.119737  0.042840  0.000000  0.745248  0.000000  0.000000  0.000000  0.000000
    0.980098  0.044852  0.000000  0.254328  0.042052  0.000000  0.000000  0.000000
    0.531336  0.056528  0.000000  0.000424  0.000322  0.044589  0.000000  0.000978
    1.058029  0.117896  0.000000  0.000000  0.000000  0.000000  0.052246  0.000000
];

numg = 4;  numm = 7;
w=21.608; % assembly width in cm
xcm  = [  0  0.5*w  1.5*w  2.5*w  3.5*w  4.5*w  5.5*w  6.5*w  7.5*w  8.5*w];
xfm  = [   10     20     20     20     20     20     20     20     20     ]/1;
ycm  = [  0  0.5*w  1.5*w  2.5*w  3.5*w  4.5*w  5.5*w  6.5*w  7.5*w  8.5*w];
yfm  = [   10     20     20     20     20     20     20     20     20     ]/1;
      
mt   = [  1  3  1  3  1  2  1  4  7  
          3  1  3  1  2  1  6  4  7
          1  3  1  2  1  3  4  7  7
          3  1  2  1  3  5  4  7  7 
          1  2  1  3  1  4  7  7  7
          2  1  3  5  4  7  7  7  7
          1  6  4  4  7  7  7  7  7
          4  4  7  7  7  7  7  7  7
          7  7  7  7  7  7  7  7  7 ];      
BCL = 1;
BCR = 0;
BCB = 1;
BCT = 0;        

src(:,:,1) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ];            
src(:,:,2) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ];           
src(:,:,3) = [ 0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0
               0  0  0  0  0 
               0  0  0  0  0 ]; 

IBSL = [0 0 0 0 0
        0 0 0 0 0]';
IBSR = [0 0 0 0 0
        0 0 0 0 0]';
IBSB = [0 0 0 0 0
        0 0 0 0 0]';
IBST = [0 0 0 0 0
        0 0 0 0 0]';